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I. INTRODUCTION 



II. GENERAL SCALAR-TENSOR THEORY 
AND ITS NEWTONIAN LIMIT 



Modern cosmological observations, like galaxies sur- 
veys (SDSS, 2dF), galaxy rotation curves, the Bullet 
Cluster observation, studies of clusters of galaxies, sur- 
veys of supernovae la, the cosmic microwave background 
radiation, and the primordial nucleosynthesis establish 
that the Universe behaves as dominated by dark matter 
(DM) and dark energy. However, the direct evidence for 
the existence of these invisible components remains lack- 
ing. Several theories that would modify our understand- 
ing of gravity have been proposed in order to explain the 
large scale structure formation in the Universe and the 
galactic dynamics. 

During the last decades there have been several propos- 
als to explain DM, for example: Massive Compact Halo 
Objects (Machos), Weakly Interacting Massive Particles 
(WIMPs). Other models propose that there is no dark 
matter and use general relativity with an appropriate 
equation of state. Or we can use scalar fields, minimally 
or non minimally coupled to the geometry. 

In this work we are mainly concern with the study of 
dark matter (DM) and its influence on the rotation curves 
of spiral galaxies. Scalar fields are the most natural can- 
didates to model dark matter models. Our DM model is 
based on using a scalar field (SF) that is coupled non- 
minimally with the Ricci scalar in the Lagrangian that 
gives Einstein field equations. 

On galactic scales we can test dark matter models us- 
ing observed rotation curves. Nowadays we have lots of 
compiled samples of galaxy data. The best type of galax- 
ies to test dark matter model is the low surface brightness 
type of galaxies. In this work we test our dark matter 
model using a sample of low and high surface brightness 
galaxies. 

We organize our work in the following form: In the 
next section we present the general theory of a typical 
scalar-tensor theory (STT), i. e., a theory that general- 
izes Einstein's general relativity by including the contri- 
bution of a scalar field that couples non-mininially to the 
metric. Next we present our results for the estimation 
of the parameters of the scalar field dark matter model. 
Finally, our conclusions are given in the last section. 



We start with the Lagrangian of a general scalar-tensor 
theory 
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Here g^j^ is the metric, Ciii{g^v) is the matter Lagrangian 
and w(0) and l/(<^) are arbitrary functions of the scalar 
field. The fact that we have a potential term V{(j)) tells 
us that we are dealing with a massive scalar field. Also, 
the first term in the brackets, is the one that gives 
the name of non-minimally coupled scalar field. 

When we make the variations of the action, S — 
J d^x £, with respect to the metric and the scalar field 
we obtain the Einstein field equations [T] 
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for the metric g^t^ and for the massive SF (p we have 



□ 0: 



1 



3 + 2a; 



[SnT - uj' {d(t>f + (l)V' - 2V] 



(3) 



where ()' = g^. Here T/^^i, is the energy-momentum ten- 
sor with trace T, uj{(j)) and V{(f>) are in general arbitrary 
functions that govern kinetic and potential contribution 
of the SF. The potential contribution, V{cl)), provides 
mass to the SF, denoted here by msp- 



A. Newtonian limit of a STT 

To study the rotation curves of galaxies we need to 
consider the influence of SF in the limit of a static STT, 
and then we need to describe the theory in its Newto- 
nian approximation, that is, where gravity and the SF 
are weak (and time independent) and velocities of dark 
matter particles are non-relativistic. We expect to have 
small deviations of the SF around the background field, 
defined here as (0) and can be understood as the scalar 
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field beyond all matter. Accordingly we assume that the 
SF oscillates around the constant background field 
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where 77^1/ is the Minkowski metric. Then, Newtonian 
approximation gives [2HS] 
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we have set {(j)) = (1 + a)/GM and a = 1/(3 + 2w). In 
the above expansion we have set the cosmological con- 
stant term equal to zero, since on small galactic scales 
its influence should be negligible. 

Note that equation Q can be cast as a Poisson equa- 
tion for V = (l/2)(/ioo + 0/(0)), 
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and the New Newtonian potential is given by $7v = 
(l/2)/ioo = ip — (1/2)0/(0). Above equation together 
with 
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form a Poisson-Helmholtz equation and gives 
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which represents the Newtonian limit of the STT with 
arbitrary potential V{(j)) and function a;(0) that where 
Taylor expanded around (0). The resulting equations 
are then distinguished by the constants Gat, a, and A = 
hp/mspc- Here hp is Planck's constant. 

The next step is to find solutions for this new Newto- 
nian potential given a density profile, that is, to find the 
so-called potential-density pairs. General solutions to 
Eqs. (|6| and ^ can be found in terms of the correspond- 
ing Green functions, and the new Newtonian potential is 
011] 
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The first term of Eq. ([8]), is the contribution of the usual 
Newtonian gravitation (without SF), while information 
about the SF is contained in the second term, that is, 
arising from the influence function determined by the 
modified Helmholtz Green function, where the coupling 
Lo (a) enters as part of a source factor. 



The potential of a single particle of mass m can be 
easily obtained from (|8]) and is given by 
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For local scales, r ^ A, deviations from the Newtonian 
theory are exponentially suppressed, and for r ^ A the 
Newtonian constant diminishes (augments) to Gjv/(1 + 
a) for positive (negative) a. This means that equation 
([9]) fulfills all local tests of the Newtonian dynamics, and 
it is only constrained by experiments or tests on scales 
larger than -or of the order of- A, which in our case is of 
the order of galactic scales. In contrast, the potential in 
the form of equation <i>Ar = -G™ (l -I- ae"''/^), with the 
gravitational constant defined as usual, does not fulfills 
the local tests of the Newtonian dynamics [7] . 



B. Multipole expansion of the Poisson-Helmholtz 
equations 

The solutions for a spherically symmetric distribution 
of mass is as follows: The Poisson's Green function can be 
expanded in terms of the spherical harmonics, Yin{6, Lp), 

^-4.5:5:^r,;(.',,')ii«(^,.), 

where r< is the smaller of |r| and jr^l, and r> is the 
larger of |r| and jr^l and it allows us that the standard 
gravitational potential due to a distribution of mass p(r), 
without considering the boundary condition, can be writ- 
ten as [Sj 

where i/)^*^ (V'*''^'') are the internal (external) multipole 
expansion of '0, 



1=0 n=-l 



1=0 n=~l 



Here, the coefficients of the expansions ■0^'"' and 
known as internal and external multipoles, respectively, 
are given by 

f dr'-^Y;j9\^')p{r'), 

Jv{r<r') ^ 

<lln=^[ dr'Y;^i9',ip')r"p{r'). 

JV(r>r') 

The integrals are done in a region V where r < r' for 
the internal multipoles and in a region V where r > r' 
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for the external multipoles. They have the property 
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We may write expansions above in cartesian coordi- 
nates up to the monopole. For the internal multipole 
expansion we have 
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and its force is 



where 
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For the external multipoles we have 



and its force is 
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The external monopole have the usual meaning, i.e., 
M^^) is the mass of the volume V{r > r'). We may atach 
to the internal monopole similar meaning, i.e., M^'^ is 
the internal "mass" of the volume V{r < r'). 

In the case of the scalar field, with the expansion 
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the contribution of the scalar field to the Newtonian grav- 
itational potential can be written as 
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where, for simplicity of notation, we are using m 
msF = hp/{c\) and 
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ii{x) and ki{x) are the modified spherical Bessel func- 
tions. 
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We have defined the multipoles for the scalar field as 
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They, also, have the property 
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The above expansions of SF contribution to the New- 
tonian potential can be written in cartesian coordinates. 
The internal multipole expansion of the SF contribution, 
up to the monopole is 
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and its force is 
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In the exterior region the SF monopole contribution to 
the potential is 
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and its force is 
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In the limit when m — >■ we recover the standard New- 
tonian potential and force expressions. 



III. RESULTS 

In this section we will show how to obtain values of 
parameters of the model, i.e., a and A using observed 
rotation curves of galaxies. 

Our SF galaxy model will be as follows. A test particle 
will move under the action of the potential 

(/){r) = (j>D (r) + (l)NDM (r) + (t>sF (r) (24) 

where (I)d{t) is the potential due to the disk mass density 

profiles of gas and stars, p^,. 4>NDM{f) is the Newtonian 
potential due to a DM mass distribution, and ipspir) is 
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the potential due to the SF contribution. They 
the equations: 


satisfy 


Galaxy 
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Sm 




1.6 


0.8846 






DDO 47 


IBm 


0.1 
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0.01 
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NGC 3109 


SBsm 




1.2 


0.11 
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ESQ 116-G12 SBm 


4.6 


1.7 


2.1 




NGC 7339 


SABbc 


7.3 


1.5 
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M31 


SAb 


0.1 

20 


0.5 
4.5 


0.01 
126 


The rotational velocity is obtained using the 
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sion V = r\d4>/dr\. For an exponential disk we have the 
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UGC 10981 
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5.4 
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where loiy), Kq^jj), and Ki{y), are the modified 

Bessel functions. Mo is the total mass of the disk and 
Ru is its scale length. 

We consider that the DM mass density profile is given 
by the pseudo isothermic profile, 
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where Ps is the central density of the halo and is its 
radius. This density profile has the force contribution 
given by 
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with the auxiliary functions, 
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TABLE I. Properties and parameters of the analyzed sample. 
From left to right, the columns read: name of the galaxy; 
Hubble clasification; Luminosity; disk scale length in kpc; disk 
mass m 



and the rotation curve is obtained from 
V\r) = VUr) + VlAr) 
We best fit the observations with 
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Free parameter are: r^, Ps, a, and A. Units we are using 
are such that a — 1 kpc = 1, Vb = 1 km/s = 1 and G = 1. 

In table |T] we show properties and parameters of the 
analyzed sample of galaxies. This sample include low 
and high surface brightness: UGC 4325 is a late-type 
dwarf galaxy at a distance of 10.1 Mpc; has an abso- 
lute B-band magnitude of Mb = -17.2 [10]. DDO 47 
is a nearby dwarf galaxy of type IBm [TT]. NGC 3109 
is a nearby Magellanic-type galaxy with a Mb — —16.35 
|12) . ESO 116-G12 is a SBm type galaxy with magnitud 
Mi = -20.0 Galaxy NGC 7339 is a SABbc type 

galaxy with magnitud Mj = -20.6 [13]. NGC 6822 is a 
nearby Local Group member galaxy with a _B-magnitude 
of Mb = —15.8 [12]. Andromeda galaxy M31 is a very 
large local group spiral galaxy, SAb type, at a distance of 
0.78 Mpc [H]. UGC 8017 is a Sab type galaxy at a dis- 
tance of 102.7 Mpc [11]. UGC 11455 is a Sc type galaxy 
at a distance of 75.4 Mpc [15]. UGC 10981 is a Sbc type 
galaxy at a distance of 155 Mpc [T?j. 

We have divided the sample in two groups. Group 
A is composed of galaxies: UGC 4325, DDO 47, NGC 
3109, ESO 116-G12, NGC 7339. Group B is composed 
of galaxies: NGC 6822, M31, UGC 8017, UGC 11455 and 
UGC 10981. 

In table |ll] we show the fitting results. We observe 
that for the galaxies in group A (UGC 4325, DDO 47, 
NGC 3109, ESO 116G12 and NGC 7339) the parameter 
a of the scalar field dark matter model is negative and 
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FIG. 1. Group A of analyzed galaxies: UGC 4325, DDO 47, piQ. 2. Group B of analyzed galaxies: NGC 6822, M31, UGC 
NGC 3109, ESO 116-G12, NGC 7339. 8O17, UGC 11455, UGC 10981. 
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ps (po) 


ra (kpc) 


A (kpc) 
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Group A: 












UGC 4325 


3492.04 


1.95 


0.691 


-0.477 


7.20 


DDO 47 


1281.99 


2.983 


1.976 


-0.259 


3.35 


NGC 3109 


800.66 


2.98 
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-0.26 
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NGC 7339 
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NGC 6822 


3501.96 


1.46 


1.09 


0.67 


2.21 


M31 


40001 


0.98 


1.79 


1.05 


3.50 


UGC 8017 


26000 
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3.09 
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UGC 11455 


4000.7 


6.30 
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0.42 
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35001 
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TABLE II. Properties and parameters of the analyzed sample. 
From left to right, the columns read: name of the galaxy; 
central density in units of po = 18491.4 solar masses per kpc'^; 
central radius in kpc; scalar field scale length in kpc; scalar 
field strength; the Xred value. 

a = —0.303 on the average. For this group the average 
value of scale length of the scalar field is A = 1.72 kpc. 
For all the rest of galaxies, group B, a is positive with an 
average value a — 0.934. In this case all galaxies, except 
NGC 6822, have a high luminosity. For this group the 



average value of scale length of the scalar field is A = 1.81 
kpc. 



IV. CONCLUSIONS 

We have used a general, static STT, that is compatible 
with local observations by the appropriate definition of 
the background field constant, i. e., (0) = (1 + a)/GN, 
to study rotation curves of spiral galaxies. 

It is important to note that particles in our model are 
gravitating particles and that the SF acts as a mecha- 
nism that modifies gravity. The effective mass of the SF 
{msF = 1/A) only sets an interaction length scale for the 
Yukawa term. 

We have estimated the parameters of the scalar field 
dark matter model by minimizing the appropriate for 
two samples of observed rotation curves. For the galaxy 
group A, i.e., the dwarf and low surface brightness galax- 
ies, the strength parameter has a negative value, and on 
the average a — —0.303 and A = 1.72 kpc. Whereas, 
for group B, i.e., the high surface brightness galaxies, 
a = 0.934 and A = 1.81 kpc. We have found, in general, 
also higher values of Xred group B of galaxies. 

This should be consistent with the general believe that 
low surface brightness galaxies are dominated by dark 
matter whereas high surface brightness galaxies are not. 
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